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ABSTRACT 

Recently Disney et al. (2008) found a striking correlation among the five basic parameters that govern the 
galactic dynamics: R^, Rgo, L r , Mj, and Mhi- They suggested that this is in conflict with the ACDM model, 
which predicts the hierarchical formation of cosmic structures from bottom up. Considering the importance of 
the issue, we performed a similar analysis with a significantly larger database and two additional parameters, 
Lj and Rj, of the near-infrared J band. We used databases from the Arecibo Legacy Fast Arecibo L-band Feed 
Array Survey for the atomic gas properties, the Sloan Digital Sky Survey for the optical properties, and the 
Two Micron All Sky Survey for the near-infrared properties, of the galaxies. We found that the first principal 
component dominates the correlations among the five parameters and can explain 83% of the variation in the 
data. When color (g-i) is included, the first component still dominates and the color forms a second PC that is 
almost independent of other parameters. The overall trend in our near-infrared PCA is very similar, except that 
color (i—J) seems even more decoupled from all other parameters. The dominance of the first PC suggests that 
the structure of galaxies is governed by a single physical parameter. This confirms the results in Disney et al. 
However, based on the importance of the baryon physics in galaxy formation, we find it premature to conclude 
that the hierarchical structure formation scenario and the notion of cold dark matter are necessarily flawed. 
Subject headings: galaxies: fundamental parameters — galaxies: structure — galaxies: statistics — galaxies: 
formation — cosmology: dark matter 



1. INTRODUCTION 

One way to understand our universe is to gain insights into 
the structure of galaxies. For one thing, it helps to reveal the 
role of dark matter in their formation and dynamics. The cos- 
mological model that consists of a cosmological constant and 
the cold dark matter in addition to the ordinary baryon mat- 
ter and radiation (ACDM) has been able to successfully ex- 
plain the evolution of the cosmic structures espec ially at large 
scales . By measuring CMB fluctuatio ns (COBE:ISmoot et alj 
1991 iBeniiett et all 1 1 9941) ( WMAP: IBennett et all l2003allbl: 
Kmatsuo et alj 120091 iDunklev et ail 120091) , Type l a super- 
novae (Riess et a 1J I1998D and gravitational lensing (IZwickyl 
119371) . this model of cosmology has by now been established 
as the standard model of cosmology. Among its various impli- 
cations, it suggests a hierarchical, bottom-up history of struc- 
ture formation that evolves from small fluctuations to galax- 
ies, clusters and eventually superclusters. 

At small scales, the success of ACDM has not been as 
clearcut. There still exist several inconsistencies between 
ACDM and the observations at small scales. For example, 
the simulations based on ACDM have revealed more number 
of galactic satellites (fMoore et al. l ll999t iKlvpin et all 119991: 
lPrimackl2009l lKravtsovl20 1 Oh and less number of disk galax- 
ies dSommer-Larsen et alj 12001b than what have been ob- 
served. Besides, the degree of emptiness in t he voids is als o 
inconsistent between theory and observation (Peeblesl l2001l) . 
While the ACDM can exp lain th e galaxy rot ational curves at 
large radii dde Bloket al.lll997t iFriedmannl [20091) . the rela- 
tively higher density at the galactic core than that predicted 
by the ACD M, that is, the so-called cusp-core problem, is still 
unresolved dDalcanton et al.ll2001l) . Aside from these, there 
is the newly emerged challenge to the CDM model with re- 



gard to the hierarchical galactic structure formation. That is, 
if galaxies are the results of random hierarchical merger pro- 
cesses, their basic parameters are expected to be independent. 
Yet a recent study found that the basic parameters of galaxies 
are highly correlated. 

The hierarchical galaxy formation scenario ha s been 
actively investigated previou sly (e.g., IWhite etafl I199U 
ICole et al1l2000t lBaughll2006[) . and the correlation between 
selected galactic parameters has bee n studied in the past 
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20081) . However, the overall correlation be- 



tween all major ga lactic variables has n ever been investigated 
until recently by dDisnev et al.1 120081 hereafter D08), who 
made the remarkable finding, based on 195 galaxies, that there 
exist very strong correlations between six galactic variables: 
the 90%-light radius (fl 90 ), the 50%-light radius (R 50 ), the H I 
mass {Mhi), the dynamical mass (M^), the luminosity (L), and 
the color. They suggested that this was in conflict with the 
notion of the hierarchical structure formation scenario, which 
would predict more collisions among substructures and there- 
fore less regularities between galaxies. Since CDM model is 
the underlying tenet of the hierarchical scenario, the authors 
assert that CDM model, or even the notion of dark matter it- 
self, is in doubt. 

Considering the importance of this issue, we set out to rein- 
vestigate the analysis made by D08. We note that their claim 
actually involves two separate sub-issues. First, there is the 
issue of whether the galactic dynamics is indeed controlled 
by a single parameter. Second, even if this is true, there still 
remains the issue of whether such a fact necessarily concludes 
the failure of the hierarchical structure formation scenario and 
the cold dark matter model. In our view the positive conclu- 
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sion of the former does not necessarily lead to the affirmative 
conclusion of the latter. In this paper we demonstrate tentative 
evidence that there may exist more than one principal compo- 
nent among the global parameters of galaxies with regard to 
the first issue. As for the second issue, we argue about the 
importance of the non-gravitational baryon physics in galac- 
tic structure formation, which renders the naive extrapolation 
of the hierarchical structure formation scenario from cosmic 
to galactic scales questionable. 

In order to scrutinize the correction issue, we performed a 
similar analysis on global parameters of galaxies with a sig- 
nificantly larger database and two additional parameters, Lj 
and Mj, based on the infrared J band. We include a total of 
1022 galaxies from the Arecibo Legacy Fast Arecibo L-band 
Feed Array Survey (ALFALFA; [Stavelev-Smith et alj 1 1996b 
IGiovanelli et alj2007tlSaintonge et al.l2008tlKent et all2008l) 
for the ato mic gas properties, and t he Sloan Digital Sky Sur- 
vey ( SPSS; lAbazauan etail 120091: lAdelman-McCarthv et all 
120091) for the optical properties. Among the 1022 galaxies, 
48 1 of them have also been detected by the Two Micron All 
Sky Survey (2MASS; ISkrutskie etail 120061) . We use these 
galaxies to study their near-infrared properties. 

Recently, several large blind 21cm surveys for galaxies 
have been conducted. Extragalactic H I surveys provide a 
new way of investigating galaxy properties other than the 
stellar luminosity, surface brightness and morphology. Be- 
cause of the largeness of the database, the co rrelations be- 
tween H I and other propert ies can be studied dHavnes et alj 
1 19841: [Rosenberg et al.l2005l) . The H I observation has by now 
been established as a useful tracer to study galaxies, especially 
those with low surface brightness. In addition, the dynamics 
of the neutral hydrogen gas in principle constrains the dark 
matter content of the galaxies. It therefore makes sense to in- 
clude the H I data in the galactic parameter correlation studies 
with the hope that it would provide us additional information 
about galaxy formation and evolution. 

Near-infrared studies of H I sele cted galaxies had been 
attempted by various grou ps (e.g., iRosenberg et al.l 120051 : 
iGarcia-Appadoo et alj|2009r) . but with relatively small sam- 
ples. Our motivation of adding the near-infrared data in our 
analysis is the following. The optical emission is sensitive to 
young stars. The near-infrared emission, on the other hand, is 
less affected by the young stars and is therefore a better tracer 
of the total stellar mass, which dominates the baryonic matter 
at galactic scales. We believe that the inclusion of the infrared 
data would provide us additional and independent information 
on the baryonic mass assembly history of the galaxies. 

Thr ough the principal component analysis (PCA; iJolliffd 
119861) . we confirm that except the color, all other observables, 
from H I, optical to near-infrared bands, are highly correlated 
and dominated by a single parameter. This is true both in the 
optical and in the near-infrared bands and this confirms the 
results in D08. However, based on the importance of the non- 
gravitational baryon physics in galaxy formation, we are not 
convinced that the hierarchical structure formation scenario 
and the notion of cold dark matter are necessarily flawed. 

The organization of this paper is as follows. We describe 
the data and sources in §2. Then several variables are adopted 
and apply to statistical analysis in §3. In §4, we summarize 
and discuss our results. 

2. DATA 

Our samples are the blind 21 cm survey from ALFALFA. 
This selects gas rich galaxies which also contain low lumi- 



nosity and low surface brightness galaxies in higher propor- 
tion than those in an optical selection. The optical data for 
this study come from the SDSS DR7, which covers 12,000 
deg 2 for imaging and provides spectra of 930,000 galaxies. 
Here we briefly describe how we selec t SDSS counterparts to 
the AL FALFA sources, and we refer to lGarcia-Appadoo et alJ 
J2009h for more detailed discussion about identification. 

The Arecibo Telescope has a beam size of 3 '5 at 21 
cm. Since the majority of the H I detections have S/N > 
10 (IGiovanelli et alj 120071: ISaintonge et alj 120081: [Kent et alj 
120081) . we adopted a conservative searching radius of 10". 
We found 1233 SDSS galaxies that appear to be detected by 
ALFALFA. We then excluded Virgo galaxies, because their 
neutral hydr ogen is known to u ndergo strong environmental 
impact (e.g., IChung et alj|2009h . We also excluded galaxies 
whose half-light radii are too small (< 1") comparing to the 
SDSS resolution. Such small half-light radii either arise from 
mis-identifications (from stars) or would result in large uncer- 
tainties. We are left with a large sample of 1022 SDSS galax- 
ies. Among these galaxies, 889 have g magnitudes that are 
< 18 and 120 have g = 18-20. The cumulat ive number counts 
of SDSS galaxies (e.g. jYasuda et al.ll2001l) are ~ 60 deg -2 at 
g < 18 and ~ 450 deg" 2 at g = 18-20. Given our 10" search 
radius, we therefore expect at most 1.3 mis-identifications in 
our 889 g< 18 galaxies, and additional 1.3 mis-identifications 
in our 120 g = 18-20 galaxies. These number are sufficiently 
small and mis-identified galaxies should not impact our anal- 
yses. 

The total stellar masses of galaxies are more directly re- 
flected by the near-infrared observations. We added the data 
from 2MASS to our samples. 2MASS provides J, H, and K s - 
band observations of the entire sky as well as a point source 
catalog and an extended source catalog. We use the 2MASS 
All-Sky Extended Source Catalog (XSC) to find the galaxies 
in the H I samples. To understand the quality of the identifica- 
tion, we first compared the 2MASS and SDSS coordinates of 
the 2MASS detected galaxies. We found that more than 90% 
of them have offsets between 2MASS and SDSS that are well 
within their half-light radii. We visually inspected all galaxies 
with large offsets of > 2", and found small number of cases 
that are likely mis-identifications as well as ongoing merg- 
ers. We excluded these galaxies from our samples. Because 
2MASS is shallower than SDSS, we are left with 481 reliably 
identified galaxies in the near-infrared. 

From the ALFALFA released catalog 1, 2 and 3, we ob- 
tained 1796 H I data, out of which 1265 galaxies could be 
found in the SDSS DR7 database. There are 32 galaxies 
within this set that are too faint in the optical to have reli- 
able magnitudes and luminosities. Hence we finally used the 
remaining 1022 galaxies in our analyses. We h ave also an- 
alyze d the 195 galaxies of D 8 from HIPASS dMeveret alJ 
2004t lAbazaiian et all I2004t IGarcia-Appadoo et al l 120*091 
West et alJ 120091 1201 Oh and have obtained similar results. 
However, since the definitions of the observational variables, 
such as that of the rotational velocity, are not entirely con- 
sistent between the two data sets, we only report on the re- 
sults from ALFALFA. These 1022 galaxies can be regarded 
as a blind H I selected sample. We deduced from the data 
six variables, which are R^ (half-light radius in units of pc) 
Rgo (90%-light radius in units of pc), L r (luminosity in r band 
in solar units), Mm (H I mass in solar units), Mj (dynamical 
mass in solar units) and color (g-i). 

The variables R50 and R90 represent the radii in the 
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Petrosian system dPetrosianl Il976t iBlanton etal] 120011: 
lYasuda et aT]|2001l) . In SDSS, the parameters are petroR5 
and petroR90, respectively. Because the Petrosian system 
is based on circular objects, we further corrected the radii 
with the major-to-minor axis ratios, which are the parame- 
ters deVAB_r or expABr in S PSS. To do this, we fol - 
low the result in IWestl d2005l) and IWest et al.l d2009l 120101) . 
The authors fitted the corrections from circular to elliptical 
apertures as functions of major- to-minor axis ratios. We di- 
rectly adopted their formulas for our corrections. By com- 
paring the likelihoods of the de Vaucouleur and the expo- 
nential models, we chose the one with the larger likelihood 
between deVAB_r and expAB_r. L r was derived from 
the Petrosian system and calculated from the Petrosian mag- 
nitude, petroMag_r. Color is from the model magni- 
tude, modelMag_g and modelMag_i, which is suggested 
by the SDSS database. The variable M4 is calculated by 
(AV) 2 Rgo/G, in which AV is the rotational velocity from the 
H I spectra and corrected for the inclination with the major- 
to-minor axis ratio as what we did for R50 and Rgo. Mhi is 
acquired directly from the ALFALFA database and it is de- 
rived from the H I flux. 

Because there is high degeneracy between J, H, and K s 
bands, we only chose J band to represent the near-infrared 
data. Therefore, we acquired R j (half-light radius in J band 
in units of pc) and Lj (luminosity in / band in solar units) 
of the 481 galaxies that overlap in the ALFALFA, SDSS and 
2MASS catalogs to gain a better insight into the stars of the 
galaxies. More specifically, Lj is derived from the magnitude 
in / band, which is the parameter j_m_ext in the 2MASS 
database. This magnitude is based on the extrapolation of the 
fit to the surface brightness profile. And Rj is the integrated 
half-flux radius of J band, which is the parameter j_r_e f f 
in the 2MASS database. 

3. METHODS AND RESULTS 

Our sampl e of 1022 galaxies is n ot onl y larger than that 
in D08 and iGarcia-Appadoo et al.l ( 120091) . but also covers 
broader ranges of size, luminosity, and mass (Figure[T|i. The 
minimum value of L r in our sample is much smaller than that 
of D08. As for the H I and the dynamical masses, although the 
median values of the two samples are similar, our sample con- 
tains a substantial amount of lower mass galaxies. Our sample 
also covers a broader range of the g — r color. Because of the 
larger sampling and the wider range of the galactic properties, 
our sample is in general more representative. However since 
2MASS is shallower, our 2MASS subsample of 481 galaxies 
is relatively speaking less representative than that from SDSS. 
Even so, our 2MASS subsample is still substantially larger 
than that of D08. 

It is important to investigate whether our H I selected sam- 
ple is biased against early type galaxies, since such galaxies 
are usually gas poor. To do so, we identify spheroidal and 
disk galaxies in our sam ple based on the morp hology with a 
method similar to that in lStrateva et all (120011) . In the SDSS 
database, there are de Vaucouleur and exponential models for 
each galaxy. By comparing the likelihood and the fractions 
of the two models for our 1022 galaxies, we found that 804 
galaxies are disk-like and 218 galaxies are spheroidal-like. 

In Figure [2] we show a color-luminosity diagram for our 
1022 galaxies. The spheroidal galaxies are in general more 
luminous, and redder than the disk galaxies. This is con- 
sistent with what we expect for elliptical and spiral galax- 
ies. Most importantly, in the color-luminosity space, the 
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FIG. 1. — Distribution of the variables. The solid histograms are our 
samples and the dashed histograms are the samples of D08. We could 
only find 157 galaxies wi th rotational velocities for Mj in D08 and 
IGarcia-Appadoo et al.1 J2009h . 

spheroidal galaxies are redder than the blue cloud although 
they do not yet form a complete red sequence. Our sample 
thus appears to include a fair number of red and elliptical 
galaxies. Although the bias against extremely gas-poor galax- 
ies can be hardly avoided here, fortunately, we found no major 
difference between these two types in our subsequent studies. 
We therefore believe that the omission of extremely gas-poor 
galaxies should not have caused major systematic bias in our 
analysis. 



5 



4 




FIG. 2. — The u — i vs. L r diagram. Gray dots are disk galaxies that are 
better described by exponential models. Black dots are spheroidal galaxies 
that are better described by de Vaucouleur models. 

Relations among the key parameters can be inferred from 
the 1022 galaxies in both ALFALFA and SDSS. For instance, 
it is found that the half lig ht radius is proport ional to the 
90%-light radius (R 50 oc fc>n: lvan der Kruitlll988l) : the r-band 
luminosity is proportional to the cubic power of the half- 
light radius (L r oc Rl ; IBlanton. et al.l 120031) : the H I mass 
is proportional to the square of the half-light radius (M HI oc 
R%: lHavnes etal] [19841: iRosenberg et al J 120051) . and finally, 
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the dyna mical mass is proport ional to the r-band luminosity 
(M d oc L r ; iGavazzi et alj[l996l) . We found that all the corre- 
lations are evident even after including the near-infrared vari- 
ables, except the color (Figure[3]and Figure|4|. 
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FIG. 3. — Scatter plots showing correlations between six measured vari- 
ables. All the variables are in solar units and with logarithmic representation. 
The diagonal line is the histograms, which have vertical scales from to 600. 
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FIG. 4. — Scatter plots showing correlations between eight measured vari- 
ables, including 2MASS data and reducing to 481 galaxies. All the variables 
are in solar units and with logarithmic representation. The diagonal line is 
the histograms, which have vertical scales from to 250. There are small 
numbers of outliers in many of the plots. They are likely mis-identifications 
or bad photometry, and they do not impact our analyses. 

As a whole, the correlations between color and other vari- 
ables are much weaker than other correlations. We tested 
various combinations of colors and found that g - i gives a 
larger PC A correlation coefficient than other colors (e.g., g— r, 
adopted by D08). The reason could be the larger wavelength 



difference between g and i. Among the three 2MASS bands, 
the result based on / is somewhat better in the PCA, pos- 
sibly because of the better signal-to-noise ratio. Hence, we 
adopted i—J for the color when we analyzed the 481 galaxies 
in 2MASS. Nevertheless, our analyses show that all the corre- 
lation coefficients are smaller than 0.7, indicating that they are 
not so highly correlated with other parameters. Intuitionally, 
more luminous galaxies tend to be redder because their colors 
are dominated by older stars. In fact, the color is more com- 
plex than any other variable because of the bias introduced by 
the very luminous young stars. 

We conducted PCA to find correlations among the vari- 
ables. PCA typically produces a series of new variables called 
the principal components, namely PCI, PC2.., etc. The corre- 
lations between these principal components and the original 
variables then reveal the general correlations between the par- 
ticular variable and the others. In our case we found that the 
first principal component, PCI, is highly correlated with the 
six observational variables. We notice that while the color 
is less correlated to PCI, possibly because of the recent star 
formation, it is much more correlated with the second prin- 
cipal component, PC2. In addition to the investigation into 
the diagram of correlations, the eigenvalues of the correlation 
matrices of the original variables give quantitative informa- 
tion of the degree of correlations. For the 1022 galaxies based 
on SDSS, the eigenvalues of PCI through PC6 are 4.52, 0.70, 
0.37, 0.20, 0.17 and 0.02 (Figure |5), respectively, where the 
maximum possible value is 6. Based on common PCA crite- 
ria, eigenvalues larger than 1 are considered significant. We 
therefore only plot PCI and PC2. Next we conducted PCA 
without color, and we found the eigenvalues of PCI through 
PC5 to be 4.15, 0.41, 0.22, 0.20, and 0.02 (Figure©, respec- 
tively, where the maximum possible value is 5. 
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FIG. 5.— The PCA results for 1022 galaxies from ALFALFA and SDSS 
with colors. Here we only show the strongest ones, PCI and PC2, because 
other principal components are not significant by PCA criterions. PCI is well 
correlated with all the variables. In the first row, the color is still correlated 
with the other five variables and PC 1 . In the second row, the rightmost plot 
shows that the color is even more strongly correlated with a new principal 
component, PC2, which is not correlated with other variables. In this plot, 
PC1= 0.421ogff 50 + 0.441ogR 90 + 0.431ogL r + 0.421ogM w/ + 0.411ogM rf + 
0.31(g-r) and PC2= -0.341og# 50 -0.241ogfi 90 +0.121ogL r -0.241ogM«/ 
+ 0.081ogM rf +0.86(g-i). 

The above observations confirm the finding of D08. All 
the observed parameters are tightly correlated with PCI. The 
eigenvalue of PCI indicates that it can explain 83% of the 
variance in the data (when color is not included). Color itself 
forms a second principal component. This might be explained 
by the fact that the optical color tends to be strongly affected 
by recent star formation activities and thus carries extra in- 
formation that is unrelated to the global formation history of 
the galaxies. We will come back to the issue of color in Sec- 
tion 4. D08 claimed that the strong dominance of PCI implies 
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FIG. 6.— The PCA results for 1022 galaxies from ALFALFA and SDSS 
without colors. Here we only show the strongest ones, PCI and PC2, be- 
cause other principal components are not significant by PCA criterions. PCI 
is well correlated with all the variables. In this plot, PC1= 0.451og/?5o + 
0.471ogR 9 o + 0A5logL r + 0.441ogM«/ + 0.431ogM rf and PC2= -0.591ogR, 
- 0.431ogi?9o + 0.28 logL, + 0.181ogM«/ + 0.601ogM rf . 

a single physical parameter to govern the structure of galax- 
ies. However, this may be a oversimplified view of galaxies 
due to limited observational parameters. To further test this, 
we extend the PCA to the near-infrared. 

We included the 2MASS /-band radius and the luminosity 
to investigate the role of stars, which dominate the baryonic 
matter in galaxies. For the 481 galaxies detected by 2MASS, 
we found the eigenvalues of PCI through PC8 to be 5.41, 
1.17, 0.52, 0.39, 0.30, 0.16, 0.03 and 0.02 (FigureO, respec- 
tively, where the maximum possible value is 8. Conducting 
PCA without color again, we found that the eigenvalues of 
PCI through PC7 are 5.35, 0.62, 0.45, 0.31, 0.19, 0.05 and 
0.02 (Figure [8]), respectively, where the maximum possible 
value is 7. 
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FIG. 8.— The PCA results for 481 galaxies from ALFALFA, SDSS, and 
2MASS without colors. Here we only show the strongest ones, PCI and 
PC2 because other principal components are not significant by PCA crite- 
rions. PCI is well correlated with all the variables, as in Fig [6] In this 
plot,PCl=0.381ogR 50 + 0.401og/?<,o + 0.491ogL r + 0.491ogR/ +0.381ogL y + 
0.361ogM m + 0.341ogM rf and PC2= -0.591og# 50 -0.391ogR 90 + 0.361ogL, 
-0.221ogRy + 0.461ogL/ + 0.141ogM« ; + 0.291ogM rf . 

hand, in the near-infrared PCA in Figure |7J the i-J color al- 
most does not involve in the first principal component and 
itself forms a second component that is almost independent 
of other parameters. A potential issue here is the combina- 
tion of SDSS i and 2MASS J. Minor differences between the 
two photometric systems may hamper the correlations. To test 
this, we replaced the i-J color with the pure 2MASS J—K 
color, and we found identical results. In addition, the me- 
dian errors in g and i for the 1022 SDSS galaxies are 0.0051 
and 0.0055, respectively, and the median errors in i and J for 
the 481 2MASS galaxies are 0.0034 and 0.064, respectively. 
These translate to typical color errors of 0.0075 in g — i and 
0.064 in i-J. Both values are significantly smaller than the 
color dynamical ranges shown in Figures [3] and H] meaning 
that the distribution of the data is not dominated by measure- 
ment errors. We therefore believe that the lack of correlation 
between the near-infrared color and other parameters is real. 
We will discuss more on this result in Section 4. 

Finally, we feel obliged to check the Tully-Fisher relation 
because the rotational velocity and the luminosity are what we 
acquired directly. The results are shown in Figure [9] We see 
that the luminosity is proportional to the rotational velocity 
with the power-law index between 3 and 4 in both r band and 
J band, which agr ees with other results in the literature (e.g., 
iPierce et al.lfl992h . 



FIG. 7.— The PCA results for 481 galaxies from ALFALFA, SDSS, and 
2MASS with colors. Here we only show the strongest ones, PCI and PC2 
because other principal components are not significant by PCA criterions. 
PCI is well correlated with all the variables. The color is correlated with 
other variables and PC 1 in the first row as well as strongly correlated with 
PC2 as in Fig|5] In this plot, PC1= 0.371ogK 50 + 0.391ogR 90 + 0.391ogL,- + 
0.391og#, + 0.381ogL, + 0.361ogM/„ + 0.341ogM rf + 0.12((-7) and PC2= - 
0.331ogfi 50 - 0.251ogfi 90 + 0.081ogL r -O.OllogRj + 0.311ogL/ - 0.081ogM w/ 
+ O.OllogM,, + 0.85(1-7) 



The overall trends in the above near-infrared PCA are sim- 
ilar to those in the optical PCA. When color is not included, 
PCI dominates and can explain 76% of the variance in the 
data. The importance of PC2 slightly increases from 8% in 
the optical case to 9% in the near-infrared case here. When 
color is included, it forms another principal component by it- 
self. These again confirm the observations of D08 that only 
one physical parameter governs the dynamics of galaxies. 

A subtle but surprising difference between the optical and 
near-infrared PCAs is the behavior of color. In the optical 
PCA in Figure [5] although the g—i color forms a second prin- 
cipal component, it still weakly correlates with other parame- 
ters and it is part of the first principal component. On the other 
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FIG. 9. — The Tully-Fisher Relation in r band and 7 band. The scatter in 7 
band is smaller than that in r band. 



4. DISCUSSION AND CONCLUSION 

For the 195 galaxies analyzed in D08, the correlations 
among the parameters R=,q, R90, L g , Mhi and Md are obvi- 
ous. By selecting a much larger, 1022 overlapping samples 
from SDSS and ALFALFA, we also performed the PCA and 
confirmed through the high eigenvalues of the correlation ma- 
trices that the correlations are similarly strong. It follows that 
the radius, luminosity, H I mass and dynamical mass of those 
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chosen galaxies are tightly correlated. However, the color ap- 
pears to be much less correlated. This indeed complicates 
the situation and may be explained by a more sophisticated 
theory that would include, for example, the influence of re- 
cent star formation and very luminous young stars. Indeed 
such studies have already been pursued by many authors (e.g. , 
de Vaucouleurslll96U ITuUv et alj[l98 2t Istrateva et alJ l200U 
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D08 suggests that the optical color (g - r in their case) con- 
sists of two components: the "systematic" component that 
correlates with other parameters (and therefore involves in the 
first principal component), and the "rogue" component that is 
more or less random (and forms the second principal compo- 
nent). It is tempting to assume that the systematic component 
comes from the established stellar populations and is related 
to the global formation history of the galaxies, while the ran- 
dom component is related to the ongoing star formation activ- 
ities and can be short-lived events in the formation history. 
Surprisingly, our near-infrared analyses suggest a different 
story. The color we adopted in our near-infrared PCA is i—J. 
It is less affected by ongoing star formation and should reflect 
the overall stellar populations in the galaxies (and therefore 
should better correlate with other global parameters). How- 
ever, in our near-infrared PCA, this color is even more uncor- 
rected with other parameters. The reason of this is unclear to 
us, and we can only conclude that colors in this kind of studies 
are highly nontrivial. 

Putting colors aside, D08 and we found that the basic pa- 
rameters of galaxies are highly correlated and that there exists 
only one dominant principal component. This is true both 
in the optical and near-infrared bands. Based on the hier- 
archical galaxy formation assumption, the dark matter halo 
formation has been well studied through the merger tree pro- 
cess (e.g.,l| 
Cole et alT '. 



Kormendvl[T977HWh ite et al. 



199ltlLongaiil[200C 
■ iBower et alJ l2006t 

Kormendv etalT 120091: iFontanot et all l2009t iPeebles et all 
20101) . D08 believe that this scenario may not be consistent 
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with the simple relation between all basic galactic parame- 
ters because the processes of merger would break the origi- 
nal galactic structures. Although our analysis has confirmed 
D08's finding, we are not convinced by their interpretation. 

It is well-known that the astrophysics at small, galactic 
scales involves additional complex physical effects, such as 
that associated with the non-gravitational baryon physics. In- 
deed, t heories and simulations such as semi -analytical mod- 
elings (Nagashimaet al. 2005; Baugh 2006) have shown the 
important influence of baryon physics in the galactic forma- 
tion processes. Moreover, remarkable progresses have been 
made with hydrodynamical simulations in resolving some 
long-standing issues. For example, the merger driven CDM 
galaxy formation is known to produce too many elliptical type 
galaxies. By adding gas components in high resolution hydro- 
dynamic cosmological simulations, iGovernato et al.l (120071) 
showed that disk galaxies can quickly form after major merg- 
ers happen. Their simulations also reproduced the Tully- 
Fisher relation and pro vided a plausible solution to the miss- 
ing sa tellite problem. ;Mashchenko, Wadsl evT & Couchmanl 
( 120071) showed that stellar feedback can potentially solve the 
well know "cusp" problem in the dark matter density profile 



whe n the simulations have sufficiently high resolution s (see 
also iCeverino & Klvpinl 120091: iGovernatoetal] HoTol) . We 
therefore consider the results of galaxy properties the man- 
ifestation of a more sophisticated version of rich baryonic 
physics under the merger tree scenario. 

In summary, recent simulations have shown that baryonic 
physics can play extremely important roles in determining 
the structure of galaxies. It remains to be seen whether the 
rapid progress in putting baryon in synergy with dark mat- 
ter can eventually reproduce the simple correlation observed 
here. We thus find it premature to conclude that the CDM 
model is necessarily flawed based on the correlation of the 
global parameters of galaxies found so far. 
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